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Abstract. Stability of stationary solutions of parabolic equations is conventionally studied by 
linear stability analysis, Lyapunov functions or lower and upper functions. We discuss here another 
approach based on differential inequalities written for the L 2 norm of the solution. This method 
is appropriate for the equations with time dependent coefficients. It yields new results and is 
applicable when the usual linearization method is not applicable. 
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1. Formulation of the problem 

Large-time behavior of solutions to differential equations has been discussed in many publications, 
see, for example, ED, [ED, lfl3l . lfT5l . First, one has to establish the global existence of the solution. 
This is done in most cases by establishing an a priori estimate which implies boundedness of the 
solutions for all times. The usual approach to Lyapunov stability of solutions is to linearize the 
problem and prove that the spectrum of the linearized operator lies strictly in the left half-plane of 
the complex plane. 

In recent papers [31, ifTOl . [fTTTl . |[T2~1 . a novel approach to the stability and long-time behavior of 
solutions to abstract differential equations is developed. This approach is applied here to nonlinear 
systems of interest in biology. 

Consider the semilinear parabolic system of equations 



"Corresponding author. E-mail: volpert@math.univ-lyonl.fr 



1 



A. Ramm, V. Volpert 



Convergence of time-dependent Turing structures 



On 

— = D(t)Au + F(u,x,t) (1.1) 

in a bounded domain f2 C R. M with a sufficiently smooth boundary and with the homogeneous 
Dirichlet or Neumann boundary condition 

Ou 

u\ dQ = or — \an = 0, (1.2) 



On 



and the initial condition 



u(x, 0) = uq(x). (1.3) 

Here u = (u\, F = (Fi, F n ), D is a diagonal matrix with positive diagonal elements 
di = di(t), which can depend on t, and 

F(0,x,t) = 0, \/x Ett,t>0. (1.4) 

The vector-function F is assumed to satisfy the estimates 

sup \F(u,x,t)\ < Mi, (1.5) 

|F(v,s,t) -F(w,a;,s)| < c F (|t-s| + (1.6) 

where c F > is a constant independent of t> , u>, x, t, and F is continuous with respect to x. Under 
these conditions, u — is a stationary solution of problem (11.11 ). (11.21 ). In the examples considered 
below we assume that Ai < 3. 

Consider the operator linearized about this solution: 



L t v = DAv + F^(0,x,t)v 
acting in the Hilbert space L 2 (f2) with the domain 

Ou 

V = {ue H\Q), u\ dn = or — | sn = 0}. 

on 

Here t is considered as a parameter. 

Suppose that the spectrum of the operator L t is located in the half-plane Re A < <j(t). If 

cr(t) < d < 0, t > 0, 

then solution of problem (11.1I) - (11.3I) converges to the stationary solution u — 0. This means that the 
stationary (equilibrium) solution is exponentially stable, i.e., the solutions with sufficiently small 
initial data converge to the stationary solution u = at an exponential rate. 

A proof of this assertion is well known in the case of the abstract evolution problem of the type 

u = Au + B(t,u), u(0) = u , (1.7) 
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where A is a linear bounded operator in a Banach space, with the spectrum that lies in the half- 
plane Rez < o" < 0, and B(t, u) is a nonlinear operator satisfying the assumption 

\\B(t,u)\\<c (t)\\ U \\ p , p>l, 

where c (t) satisfies a suitable smallness assumption (see, e.g., BU). 

In O, Theorem 1.4.1, the following result is proved: if A is a bounded operator in a Banach 
space, then there exists the limit 

In ||e A *|| 

k := lim — = max{i?eA|A G a (A)}, 

t— >co t 

where a (A) is the spectrum of A. Therefore, if all the solutions to the Cauchy problem 

u = Au, u(0) = u , 

decay exponentially fast to zero, then the spectrum of A lies in the half-plane Re z < — n, k > 0, 
and vice versa. One should have in mind that if A is a bounded linear operator in a Hilbert space 
H, such that ReA < —k, k > 0, i.e., Re(Au, u) < — n\\u\\ 2 , then the spectrum of A lies in the 
half-plane Rez < —k, but the converse of this statement is false if dim if > 1: even in two- 
dimensional Hilbert space one can give an example of A with the spectrum lying in the half-plane 
Rez < — k, < 0, for which the inequality Re(Au, u) < —k\\u\\ 2 does not hold. For instance, let 
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The spectrum of this A consists of negative eigenvalue A = — 1. The quadratic form for real- 
valued u± and «2 is Rc(Au, u) = —u\ — u 2 . + 3uiu 2 . This quadratic form is not negative-definite. 

If the spectrum of A does not lie strictly in the left half-plane of the complex plane, or a(t) — > 
as t — > oo, then the assertion about exponential rate of convergence to zero of the solutions to the 
Cauchy problem (11.71 ) is not valid, in general, and the Lyapunov stability problem cannot be solved 
by a study of the linearized problem. 

In this work we study this, more difficult, case, and use a new technical tool for such a study, 
see Lemma 2.1. Let us emphasize that a(t) will not necessarily be assumed negative in this paper 
(see 021). 



2. Convergence of solutions 

In what follows we assume that F(u, x, t) satisfies assumptions made in Section 1, see (11.5b and 
(11.6b . The initial and the boundary conditions satisfy the compatibility conditions, u (x) = on 
<9f2 for the Dirichlet and |^ = on d£l for the Neumann boundary condition. Under these (and 
some additional) conditions (see, e.g., iH) there exists a classical solution of problem (ll.lb - (11.3b . 
Let us assume that 

F(u, x, t) = A(x, t)u + B(u, x, t), 
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where 

Re (A(x,t)u,u) < -7(t)|w| 2 , VieO, i > 0, (2.1) 

and 

\B(u,x,t)\ < c {t)\u\ p , Vx G ft, t > 0, p>l. (2.2) 

Here (, ) denotes the inner product in M 3 , and \u\ 2 = Y^j=i \ u i\ 2 - We assume that the diagonal 
elements di(t) of the matrix D{t) of the diffusion coefficients satisfy the estimates 

di(t) > d(t), i = l,...,n, t>0, (2.3) 

where d(t) is a positive function. The assumptions about the behavior of d(t) for large t will be 
specified below, in the formulation of Theorems 3.1-3.3. 

Let g(t) := \\u(-,t)\\, where || ■ || denotes the L 2 (f2) norm. We will also use the space 
with the norm || • W^, and the usual Sobolev space H 2 (Q) with the norm || ■ H^™. 

Multiplying equation (11.11 ) by u and integrating, we obtain, taking into account (|2.11 l- (|2.3l ): 

gg < -d{t)\\Wuf - 7 (t)/ + Co(t) / \u\? +1 dx. (2.4) 

Jn 

In the case of the Dirichlet boundary condition, we use the Poincare inequality 



c(fi) / \u\ 2 dx< / \Vu\ 2 dx, (2.5) 
Jn Jn 

where c(f2) is a positive constant which depends on the domain. The optimal (maximal possible) 
value of c(f2) is equal to the first eigenvalue Ai of the Dirichlet Laplacian in ft. 

In the case of the Neumann boundary condition, we put c(ft) = 0. 

Using the following multiplicative inequality (see, e.g., [Q]|, p. 193): 



P||oo ^ "-H"'llL2( fi )||"||^2(Q), 



|V4 || .||3/4 

where the constant c > is independent of u, we obtain 

jf \u\^dx < g*\\u\\^ < c^\\ut^g^l\ 
From this inequality, (12.41 ) and (12.51) we obtain 

g < -(d(t)c(ft) + -y(t))g + Coity-'Mf-^g^y 4 . (2.6) 

It is known that under our assumptions the if 2 norm of the solution is bounded (see [3|, Theorem 
16.1, p. 170, and Section 4). 
Define 

a(t) := d(t)c(n) + 7 (t), a(t) := c^c^ull^ 4 , g := H + l. (2.7) 
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Then (12.61 ) can be written as 

g<-a(t)g + a{t)g\ g(t) > 0, (2.8) 

where q > 1 because p > 1. 

Assume that a{t) and a(t) > are continuous functions defined on [0, oo). 

We will use in Section 3 the following basic result from iTTOl . where more general results are 
obtained (see also lfT2l ): 

Lemma 2.1. If there exists a function fi(t) > 0, defined on [0, oo), such that 

a(t) < ^ l (t) (a(t) - f^j , t > (2.9) 

and 

ti(0)g(0) < 1, (2.10) 



then g{t) exists for all t > 0, and 



0<g(t)<^- Vt>0. (2.11) 



Note that if lim^oo fi(t) = oo, then estimate (12.1 II) implies that Hindoo g(t) = 0. The function 
a(t) in lemma 2.1 is not necessarily positive. 



3. Applications 

A relatively general class of abstract differential equations for which our method is applicable is 
described by the equations of the form 

u = A(t)u + G(t,u) + f(t), u(0)=uo, 

where A is a linear operator in a Hilbert space H, G is a nonlinear operator in H, and / is a 
given function with values in H. The following assumptions allow one to use our approach: 
Re(A(t)u,u) < -70)1 Ml 2 , < a(t,g),g := ||u(t)||, ||/(t)|| < /3(t), where the func- 

tions 7,a(t, g) and (3 satisfy some assumptions that are detailed in 01, IfTOl . 

In this section we will apply the results obtained above to reaction-diffusion system (11.11) with 
time dependent coefficients. In particular, in the case where the diffusion coefficients converge to 
zero and conventional results on stability of stationary solutions are not applicable. 
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3.1. Convergence with various rates 

Exponential rate of convergence. 

In order to make clear our method for a study of the large-time behavior of the solution to 
problem (|1.11 )- (|1.31 ), let us consider first a single equation and the Dirichlet boundary condition. 
Specifically, consider the following example: 

A(x, t) = a > 0, D(t) = d > 0, c (t) = c , 
where a , d and c are some constants, and p = 2 in (12.21) . Then j(t) = —a , and 

5 

a = d c(Q) -oq, q = - . 



If a > 0, that is, if 



then we choose 



< c(ft), (3.1) 

an 



vt 



H(t) = n e 

where /i and v are positive constants. 

Let us formulate sufficient conditions for assumptions (12.91) and (12.101) to be satisfied. If these 
assumptions are satisfied, then inequality (12.1 II) yields an exponential rate of decay of the function 
g(t), and, therefore, of the solution u(t) to zero. 

This assertion can be explained in terms of the exponential stability in the sense of Lyapunov 
of the solution u = to the problem ( |l.lt - (|1.3t . Namely, consider the problem, linearized about 
the zero solution. The principal eigenvalue of the linearized problem 

doAu + aou = Ait, u\d^ = 

becomes negative if 

X < c(fi), (3.2) 
do 

where the constant c(Q) is from the Poincare inequality (12.51) . This is condition (13.11) . 
To satisfy assumption (12.101) one may choose 

Ho = 9(0)^- 

One may assume that g(0) ^ 0, because otherwise the solution is zero by the uniqueness theorem 
that holds under our assumptions. 

To satisfy assumption ( 12.91 ) it is sufficient to assume that 

Cco(t) < s(0)-M>(t) - u], (3.3) 
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where we took into account that q— 1 = | and denoted by C the constant c p_1 1 |u| • In Section 
4 it is proved that the norm | |u| \h 2 (ci) can be estimated from above by a constant independent of u. 
If one chooses v = 0.5a = (rf Ai — a )/2, then inequality (13.31) holds provided that 

co(f) < 0.5^(0)~5C _1 e^V. (3.4) 

Thus, condition (13.41) is sufficient for the assumption (12.91 ) to be satisfied. Condition (13.41) holds for 
any fixed g(0) if co(i) is sufficiently small. It holds for a fixed co(t) if g(0) is sufficiently small. 
We have proved the following result. 

Theorem 3.1. Let the function cr(t), defined in ( 12. 71) , satisfy the inequality 

a{t) > a > 0, Vt > 0, 

where <To is « constant. Choose /x(t) = /^o exp(z/t), /io = g(0) -1 , z/ = 0.5cr . If condition ( 13.41) 
holds, then the L 2 norm of the solution u(x, t) to problem ( Ii.il) -( f773l) wzY/i ?/ze Dirichlet boundary 
condition satisfies the estimate 

\H-,t)\\ < g(0)e-°- 5aot , Vt>0. (3.5) 



The conclusion of this theorem follows from Lemma 2.1, see estimate (12.1 II ). 

The method for estimating the large time behavior of solutions to evolution problems, that was 
used in the proof of Theorem 3. 1 is easy to apply in many problems. 

The assumptions of Theorem 3.1 do not explicitly require that the spectrum of the linearized 
problem lies in the open left half-plane of the complex plane. However, the exponential rate of 
decay of the solution suggests that this is the case (see 0, p. 42, p.51). 

Convergence at a power rate. 

Let us consider problem (ll.ll )- (ll.3l ) with the Dirichlet boundary condition. Let us assume that 

dif) = 7 (0 = -T^u > /*(*) = + 1)™ (3-6) 

where d(t) is the lower bound of the diffusion coefficients, see (12.31 ). do,7o, an d A*o are some 
positive constants, k > 1 is a constant. Then inequality (12.91) takes the form: 

a(t) < , t \t + 1,^-. (ctfl)^ - ^ - , t > ,3.7, 

Let us assume that 

c(Q)d >~fo + m. (3.8) 

If the above inequality holds, then the right-hand side of (13.71 ) is positive. This inequality gives a 
condition on the function c (t), defined in equation (12.21 ) and used in the definition of a(t) in (I2.7I ). 
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If m(q — 1) < 1, then condition (13.71) implies that c (t) should converge to zero as t — > oo, if 
m(q — 1) > 1, then it c (t) may grow, as t grows, and still inequality (13.71 ) may be satisfied. 
To satisfy assumption ( 12.101 ) one may choose 

fi o = g(0)-\ (3.9) 

If (13.71) . (13.81) . and (13.91) hold, then one may apply estimate (12.1 II) and obtain the following theorem. 

Theorem 3.2. If conditions H3.6\) - ( 13.91) are satisfied, then the L 2 norm of the solution u(x,t) of 
problem (li.i|) - (li.3|) with the Dirichlet boundary condition admits the estimate 

\H-,t)\\ < g(o)(t + i)- m , vt>o. 



Boundedness of the solution. 

Consider the case when global asymptotic stability of the stationary solution may not hold. We 
wish to obtain an estimate of the solution of the evolution problem, which yields stability in the 
sense of Lyapunov. We will illustrate the method in the case of Neumann boundary condition. 

If the Neumann boundary condition holds, then, in contrast with the Dirichlet boundary condi- 
tion, one has c(f2) = in equation (12.71 ) and inequality (12.51 ), so one gets a(t) = ^(t). 

Let 

7W = y^, /*(*) = AK,(l + *) m , 

where 70 and fx are some positive constants. If 70 > m and c (t) is such that inequalities (12.91) 
and (12101) hold, i.e., 

and 

then inequality (12.1 II ) yields convergence at the rate 0((1 + t)~ m ). In this example j(t) is positive. 

We can consider the case when 7(2) is negative, but then has to be a decreasing function. 
For instance, assume that 

7W = ~?J+V' /*(*) = /*o + /*i(* + l)-", (3-10) 

where the constants 7o,/io>A*i > and v > 0. In this case, (12.1 II) yields boundedness of the 
solution for all t > 0, but the solution does not converge to zero. 
Inequality (12.91 ) takes the form: 

«(t)<(Aio + A*i(* + l) ) ( ~ + ip -(TTlF J" (3 - U) 

This inequality holds if, for example, 

u+ 1 < k. 
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and, with a(t) = Cc (t), the following inequality holds: 

C(l + ty +1 c {t) < fit 1 - 70) • (3.12) 

If (T3TT2T) holds, and 

/i(0) =/i + Ail = [^(O)]" 1 , (3.13) 
then inequality ( 12.1 II ) yields the following Theorem. 

Theorem 3.3. If conditions l\3.10\) - {[3.12\) hold, and v + 1 < k, then the L 2 norm of the solution 
u(x, t) of problem di.iD -f TOl) with the Neumann boundary condition satisfies the estimate 

IK,*) II < K)] _1 <N _1 vt>o. 



3.2. Time-dependent Turing structures 

Consider a reaction-diffusion system 



3t i cP 1 1 1 

m =dl{t) d^ + F{u > v > t) > (3 " 14) 



dv , , . d 2 v 
in the interval < x < L with the boundary conditions 



<k(t)-z-2+G{u,v,t) (3.15) 



u(0) = u(L) = 0, v(0) = v(L) = 0. (3.16) 

Reaction-diffusion systems describe various applied problems, for example, biological problems. 
These systems are often considered in the case when the coefficients and the nonlinearities do not 
depend explicitly on time. We introduce time dependence in order to describe variations of the 
environment (e.g., climate factors), or to control the system behavior. For instance, if u and v are 
some concentrations, then the coefficients of mass diffusion and the reaction rates can depend on 
the temperature that can change in time due to some external conditions, or the temperature can 
serve as a control parameter. 

Suppose that F(0, 0, t) = G(0, 0, t) — for all t > 0, that is u — v — is a stationary solution 
of problem (I3.14I) - (I3.16I) . This zero solution is also a stationary solution to the ODE system 

^ = F(u,v,t), (3.17) 

(ill 

- = G(u,v,t). (3.18) 
at 

To simplify calculations, let us assume that 
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F(u,v,t) = cj){t) F°{u,v), G(u,v,t) = (f>(t)G°(u,v), d^t) = <f>(t) d?, % = 1, 2. 

Consider first the case where 0(t) = 1. Let us choose parameters in such a way that w = v = 
is a stable solution of system ( 13.171 )- (13.181) but it is unstable as a solution of problem (|3.14| )- (|3.16l ). 
In this case, another solution, which is not homogeneous in space, can appear. This is so-called 
Turing structure, that is often discussed in relation with numerous biological applications (see, e.g., 
031, The Turing structure provides one of the possible mechanisms of pattern formation 

in biology. 

We assume that the solution u — v — of system (13.171) , (13.181 ) is stable, and that the eigenval- 
ues of the matrix 



M 

have negative real parts. Here 



a b 
c d 



a = F°(0,0), 6 = 2^(0,0), c = G°(0,0), d = G° v (0,0). 

In order to study stability of the zero solution as a stationary solution of problem (I3.14I )- (I3.16I ). 
consider the linearized system 

-^ = d° 1 — + au + bv, (3.19) 

dv 0^ v 

at = d °d^ +cu+dv < 3 - 2 °) 

with the boundary conditions (I3.16I ). 

If one looks for the solutions of this system in the form 

u — pi sin(kx)e xt , v = P2 sin(kx)e xt , 
then one obtains the following eigenvalue problem: 

M(k)p = Xp, 

where 



Af(Jfe) 



a - d\k l b 



c d- d° 2 k 2 

P— (Pi)P2)» M(0) = M. Denote its eigenvalues by Aj(fc), % = 1,2. The assumption Re Aj(0) < 0, 
i — 1,2 implies: 

a + d<0, ad-bc>0. (3.21) 

Furthermore, 
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det M(k) =ad-bc- (ad° 2 + dd\)k 2 + d^k 4 , Tr M(k) = a + d-(d° 1 + d° 2 )k 2 . 

If detM(fc) = 0, then one eigenvalue of this matrix is negative and another one equals zero. 
Hence, system ( 13.191) , ( 13.201) linearized about the solution u — v — has a zero eigenvalue. If, 
under a change of parameter, this eigenvalue crosses the origin, then a spatially inhomogeneous 
solution can bifurcate from it. 

Thus, equality det M{k) = determines the stability boundary of the solution u — v — and 
the condition of bifurcation of a spatially inhomogeneous solution. 

Let us verify that equality det M{k) = is compatible with inequalities (13 .2 1 b . If a < 
and d < 0, then detM(fc) > 0. In order to have det M(k) = 0, let us assume that one of the 
coefficients a or d is positive, but the sum a + d is negative. Let us assume, for instance, that a > 0. 
The constants b, c, d can be chosen in such a way that inequalities (13.211 ) are satisfied. For some 
a, b, c, d, k, d 2 fixed, we can increase d® so that the determinant of the matrix M{k) becomes zero. 

Thus, if (p(t) = 1, then u = v = can be a stable solution of system (13.171 ), (13.181 ), but unstable 
as a solution of problem (I3.14I )- (I3.16I ). In this case, a stationary spatial structure can emerge and 
the solution of the evolution problem can converge to it. 

If cf){t) ^ 1, then the previous considerations do not allow us to conclude whether the solution 
of problem (|3.141 )- (|3.161 ) with a given initial condition converges to a trivial solution or to a spatially 
inhomogeneous solution. 

Let us use the method developed in Section 2 in order to study the behavior of solution of the 
time dependent reaction-diffusion system. We have 

<7(*) = ^(f)(d c(fi)-7o), 

where 

d = min(d 1 , d®), 



aul + Q> + c ) u i u 2 + du\ < + ^(6 + c)J u\ + ^d + ^(6 + c))u\< ^q{u\ + ul), 

7o = max + ^(6 + c),d+ ^(6 + c) 

We obtain the following result. 

Theorem 3.4. 1. Assume that d c(tt) > 7 , 

If 4> (d c(Q,) — 70) > m, and c (t) (see (12.21 )) and a(t) (see (12.71) ) are such that condition (12 .91 ) is 
satisfied, then 
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H-,t)ll <s(o)(t + i)- m , *>o. 

2. Let d c(Q) < 70, 

Kt) =fM) + ^ , (fiQ + = #(0), 

where /i , /ii and fc are some positive constants. 
If c (t) and a(t) satisfy condition (12.91) . then 

\H;t)\\<»O l Vt>0. 

The conclusion of this theorem follows from Lemma 2.1. The first part of the theorem gives a 
sufficient condition of convergence to the trivial solution. If this condition is not satisfied, then the 
solution can possibly converge to a spatially inhomogeneous solution. In this case, the second part 
of the theorem gives an estimate of the solution. 

4. An estimate of the solution 

Lemma 4.1. Suppose that for some positive constant M\ the following estimate holds: 

\F(u,x,t)\< Mi, V« e M", 16O, t>0. (4.1) 
Then solution of problem (|l.lt - (|1.3t with the Dirichlet boundary condition satisfies the estimate 

IMIi^n) < M 2 , t > 0. (4.2) 
Proof. Each component Ui of the solution satisfies the problem 

du ■ 

—^ = d l Au i + F i (u,x,t), (4.3) 
at 

Ui\ d n = 0, Ui(x,0) = u° l (x). (4.4) 
We first obtain an estimate of the solution in the uniform norm. Let 

v(x) = — a\x\ 2 + b, 

where 

-2nadi + M 1 <0. 
The ball with the radius R = yjb/a contains the domain f2, and 

u °i{. x ) < v ( x )i x E Q. 
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Such constants a and b can be chosen for any Mi and any initial condition. 
Then v(x) is an upper solution of equation (14.31 ) and 

Ui(x,t) < v(x), x e Q, t>0. 

The functions Ui(x,t) are bounded from below as well. Thus, estimate (14.21 ) follows from the 
known estimate (see, e.g., 0) 

IMI# 2 (fi) < K (ll-^ll^cn) + , 

where K > is a constant independent on u. 

□ 

Lemma 4.2. Suppose that 

Fi(u,x,t)<0, V«i><, i£fl, t>0, (4.5) 

for some constants u*, 1 < i < n. Then solution of problem (|l.lt - (|1.3t with the Neumann 
boundary condition satisfies estimate (14.21) . 

Proof. It is sufficient to note that any constant greater than u* is an upper solution of equation 
(1431) . 

□ 
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